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Abstract 

We use the solutions of the noncommutative Wheeler-De Witt equation arising from a 
Kantowski-Sachs cosmological model to compute thermodynamic properties of the Schwarzschild 
black hole. We show that the noncommutativity in the momentum sector introduces a quadratic 
term in the potential function of the black hole minisuperspace model. This potential has a local 
minimum and thus the partition function can be computed by resorting to a saddle point eval- 
uation in the neighbourhood of the minimum. The thermodynamics of the black hole is derived 
and the corrections to the usual Hawking temperature and entropy exhibit a dependence on the 
momentum noncommutative parameter, n. Moreover, we study the t = r = singularity in the 
noncommutative regime and show that in this case the wave function of the system vanishes in 
the neighbourhood of t = r = 0. 
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I. INTRODUCTION 



Black Holes (BHs) are one of the most puzzling legacies of general relativity. Their 
existence implies that space-time is geodesically incomplete and has at least one singular 
point where physical quantities are divergent. Furthermore, semiclassical considerations 
reveal, rather surprisingly, that BHs radiate and have an associated thermodynamics. It is 
therefore natural to expect that a quantum gravity theory might allow for an underlying 
statistical mechanical explanation of the thermodynamical properties of BHs and shed some 
light on the prevalence and possible evolution of space-time singularities. 

In this respect, given that a quantum gravity theory is not within grasp, a quantum cos- 
mology approach based on the mininsuperspace approximation might bring some insight 
about the prevalence of singularities under quantum conditions. Furthermore, the evidence 
that noncommutative features are present in string theory, one of the most discussed quan- 
tum gravity proposals, suggests that some noncommutativity structure should be equally 
included into the quantum cosmology approach. This means that one should seek for a 
noncommutative version of the Wheeler-De Witt equation. This generalization has been 
carried out in the context of the Kantowski-Sachs (KS) cosmological model for the most 
general (canonical) noncommutative algebra, namely the one that includes noncommuta- 
tivity of the momentum variables as well [lj . An additional bonus of this approach is that 
the KS model can be used to study the interior region of a Schwarzschild BH. 

Thus, the aim of the present work is to use the phase-space noncommutative general- 
ization of the KS cosmological model developed in Ref. lj to examine the interior of a 
Schwarzschild BH. As we shall see this will allow for a proper treatment of the BH thermo- 
dynamics, contrary to what happens when momentum noncommutativity is absent, and 
for an analysis of the impact of the quantum and noncommutative effects on the issue of 
singularities. 

Let us start our discussion with the description of a Schwarzschild BH. General relativity 
allows for solutions where the causal structure of space-time changes at different regions 
of space-time. A well known example is the Schwarzschild BH, described by the following 
metric, 

ds 2 = -(l-^dt 2 + (l-^-y 1 dr 2 + r 2 (d9 2 + sin 2 ed^ 2 ) , (1) 

where r is the radial coordinate. The fact that the first component of the metric vanishes at 
r = 2M defines the horizon of events for the Schwarzschild BH. For r < 2M the time and 
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radial coordinates are interchanged (r <->• t) so that space-time is described by the metric 

■2M \~ x , , flM 



ds=-[ 1 dt 



l^j dr 2 + t 2 (d6 2 + sin 2 6d<p 2 ) . (2) 



* / V t 

This means that an isotropic metric turns into an anisotropic one, implying that the interior 
of a Schwarzschild BH can be described by an anisotropic cosmological space-time. Indeed, 



the metric (j2J) can be mapped to the metric of the KS cosmological model 
the Misner parametrization, can be written as 



3], which, in 



ds i = -jv 2 ^* 2 + e 2 ^dr 2 + e- 2 ^e- 2 ^ n (d9 2 + sin 2 9d<p 2 ) , (3) 

where fl and f3 are scale factors, and N is the lapse function. The following identification 
for t < 2M, 

n2 = /2M _ \ 1 ^ ^ = {™L_{\ } e -2Vsi3 e -2Vsn = t 2 ; (4) 

allows for turning the metric Eq. (j3J) into the metric Eq. (T5]) away from the horizon 
t = r = 2M, the coordinate singularity. 

In what follows we shall assume that, at the quantum level, the interior of the BH is 
described by the phase-space noncommutative extension of the quantum KS cosmologi- 
cal model. In this setting, the noncommutative Wheeler-De Witt (NCWDW) equation 
has been recently obtained and its solutions calculated explicitly 1]. As we shall see, 
the introduction of momentum noncommutativity changes the WDW equation and has a 
significant impact on the functional form of its solutions. The new WDW equation al- 
lows for consistently calculating the partition function for the Schwarzschild BH through 
the Feynman-Hibbs procedure [4( and then on all the relevant thermodynamic quantities. 
This, as we will see, is in contrast with what happens in the commutative and configuration 
space noncommutative cases. Moreover, the solutions of the NCWDW equation provide 
a new insight into the problem of the BH singularity. We will show that these solutions 
vanish in the neighbourhood of t — r — 0, but that this does not necessarily imply a van- 
ishing probability of finding the system at the singularity. We conclude that the canonical 
phase-space noncommutativity discussed in this work is insufficient to ensure the avoidance 
of singularities, but quite interestingly, our results strongly suggest that the singularity may 
be removed through other, slightly different, forms of momentum noncommutativity. 

This paper is organized as follows: We start by recalling the main features of phase- 
space noncommutative quantum mechanics (NCQM) and quantum cosmology in section 



II. In section III, we use the solutions of the phase-space NCWDW equation to obtain the 
relevant thermodynamical properties of the Schwarzschild BH through the Feynman-Hibbs 
procedure. In section IV we discuss the implications of our results for the BH singularity 
issue. Section V contains our conclusions. 



II. NONCOMMUTATIVE PHASE-SPACE QUANTUM COSMOLOGY 

It is believed that the theory of quantum gravity will determine the ultimate structure 
of space-time at the Planck scale. It has been pointed out that at this scale space-time 
might have a noncommutative structure, as suggested, for instance by string theory 5|, ]6j]. 
In this work we use units such that h = c = k = G=l, which imply for the Planck length 
and mass lp = Mp = 1. 

At the quantum mechanical level, noncommutativity has been extensively discussed 
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121 ] . In this instance, one considers canonical extensions of the Heisenberg- 



Weyl algebra, where time is assumed to be a commutative parameter and the theory is set 
in a 2<i-dimensional phase-space of operators with noncommuting position and momentum 
variables. The noncommutative algebra reads, 



[qiAj] = i0 ij > [QuPj] = > \PhPj] = ir lij , i,j= l,-,d 



(5) 



where rjij and % are antisymmetric real constant (d x d) matrices and 5ij is the identity 
matrix. The extended algebra is related to the standard Heisenberg-Weyl algebra: 



Ri 1 Rj 







Rii Ilj 



ihS,, 



IJ , 



, i,j = 1, ...,d 



(6) 



by a class of linear (non-canonical) transformations: 



q { = q { (Rj, flj) , pi=pi {Rj, II, 



(7) 



In the mathematics literature this mapping is usually referred to as Darboux (D) map 111 ]. 

A phase-space noncommutative extension of a quantum cosmology was studied in Ref. 
[1]. The presence of at least two scale factors is necessary to consider such an extension 
and the KS cosmological model, based on the metric Eq. ([3]), is a most natural candidate. 
In this case, we impose the following (noncommutative) classical Poisson algebra on the 
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scale factors /?, f2 and their conjugate momenta Pp, Pq 

{Q,Pn} = l {P,Pp} = l 

(8) 

{n,P} = 9 {p n ,Pp} = v . 

The remaining brackets vanish. 

Following the ADM procedure and taking Q, (3 as configuration variables, one can derive 
the Hamiltonian constraint for this system, 



H = NH = Ne 



P 2 P 2 

f_n _|_ jj? _ 2 e -zV3n 
24 24 



(9) 



The lapse function for the Hamiltonian problem associated to metric Eq. (J3J), N, is chosen 
as N = 24e-^- 2 ^ fl. With this choice, the classical equations of motion simplify 
drastically. For the algebra ©, these read: 

(1 = -2P n , (a) 

P n = 2 v Pp - 96V3e- 2 ^ n , (b) 

l3 = 2Pp- 96VS0e- 2V * n , (c) 

Pp = 2 V P n . (d) (10) 

An analytical solution of this system is beyond reach, given the entanglement among the 
four variables. However, a numerical solution can be obtained and used to provide predic- 
tions for several relevant physical quantities Despite the difficulty to obtain analytic 
solutions one notices that Eqs. (TTOh) and (fTtH) yield a constant of motion 

Pp = -7?(-2P n ) = -r)£l Pp + rjtt = C , (11) 

which will play an important role in solving the phase-space NCWDW equation. 

The canonical quantization of the classical Hamiltonian constraint, 7i ~ 0, based on the 
ordinary Heisenberg-Weyl algebra, yields the commutative WDW equation for the wave 
function of the universe. For the simplest factor ordering of operators this equation reads 



.P 2 + p2 _ A8e -2Vsn 



tl>(n,P) = 0. (12) 



where P n = -zjj, Pp = are the fundamental momentum operators conjugate to 

Q = Q and (3 = j3, respectively. Here, on the other hand, we shall perform the quantization 
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of the classical algebra (jSj). We thus obtain the following noncommutative Heisenberg-Weyl 
algebra: 



i9 



Pn, P/3 



(13) 



IT] 



The remaining commutators vanish. 

The simplest way to implement algebra ([TBI in this context is through a suitable D map 
(i.e. a non-unitary linear transformation mapping the noncommutative algebra into the 
standard Heisenberg-Weyl algebra) []| 



Pn = vPn c + 



V 



Pc , P/3 



(14) 



2/x ' r Mc 2/i 

where from now on the index c denotes commutative variables, i.e. variables for which 

i. The transformation Eqs. (CO 







Pn c , Pp c 


= and 






"A, *k 



be inverted, provided 

£ = ^ < 1 . 

The dimensionless constants A and /i are related by [l|, 



can 



(15) 



(A/i) -A/i + - = 0^A/i = 



(16) 

The transformation Eqs. ( fl4l) provide a representation of the operators ( fTBl as self-adjoint 
operators acting on the Hilbert space L 2 (IR 2 ). In this representation the WDW Eq. (fT21) 
is deformed into a modified second order partial differential equation, which exhibits an 
explicit dependence on the noncommutative parameters 

9 d ' 



( ■ 9 ii n 



^wri^v ~ 48exp 



From Eq. (fill) we define a new constant operator A 



V 



-2VS(*l. + «- 2A8A 



(17) 



2=, from which follows that 



(18) 



This new operator commutes with the noncommutative Hamiltonian of Eq. (TTTj) . Hence 
we may look for solutions of Eq. (fT7j) that are simultaneous eigenstates of A. A generic 
eigenstate of A with eigenvalue a G 1R, satisfies 



9 77 \ 

■i^^j + t^-^c ^ a (^c,/5 c ) = aip a (£l c ,p c ) . 



(19) 



which yields 



% I Tj \ 

-la - —Q c p c 
u\ 2/i / 



(20) 



Substituting this solution into Eq. ffTTl) implies that (fi a (z) = $t, a (Q c (z)) satisfies 

pi a 

(f)" a {z) + (vz - a) 2 (t) a (z) - 48 exp [-2V3z + -\— a]<j) a (z) = , 

where we have introduced the new variable, 

Vt r d d 



(21) 



/i dz ^ dVt c 



(22) 



This is a second order ordinary differential equation that can be solved numerically. The 
equation itself depends on the eigenvalue a and on the noncommutative parameters 6 and 
V- 

As shall be seen, Eq. ( 1211) . can be used to obtain the temperature and the entropy of 
the Schwarzschild BH. 



III. THERMODYNAMICS OF PHASE-SPACE NONCOMMUTATIVE BLACK 
HOLE 

We employ now the NCWDW equation to study the quantum behaviour of the inte- 
rior of the Scharwzschild BH. Our aim is to obtain its temperature and entropy for the 
phase-space noncommutative extension. In order to get these quantities we compute the 
partition function through the Feynman-Hibbs procedure [3]. This method is based on 
the minisuperspace potential function, which can be determined from Eq. (l2Tj) : 

\/ T \f) 

V(z) = 48 exp [— 2v / 3z H — - — a }-(r]z- a) 2 . (23) 
For convenience, let us introduce a new variable, 

x = z ~ik a - (24) 

Thus, the potential function becomes, 

V(x) = 48 exp (-2y/3x) - (r]x - cf , (25) 

where c = Pp{Q) +r]Q(0) is a constant from the classical constraint. This potential function, 
which is depicted in Fig. [T]b, exhibits a local minimum, xq, and a local maximum. 
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(a) j] = and c = 0.01 (b) r) = 1.5 and c = 5.68 

FIG. 1: Potential function for some typical values of r/ and c: (a) The potential function for the 
noncommutative case, 9^0 and 77 = and (b) The potential function with 77 7^ 0. 



In Fig. [T] we present the potential function for the noncommutative case, 9 7^ and 
77 = 0, and the noncommutative case, where 77 = 1.5. The case with noncommutativity 
only in the configuration variables (i.e. 77 = 0, 9 7^ 0) [14( is thus qualitatively similar 
to the commutative one, i.e. it has no local minimum. As can be seen in Eq. ( 12 3 p . the 
noncommutative parameter 9 appears in the argument of the exponential. Thus, unless 
noncommutativity in the momentum sector is present one does not have a local minimum 
of the potential for finite x. The values of 77 used to plot the potential are fairly typical 
and the qualitative behaviour of the potential function is identical for other values. The 
constant c is obtained using Pp(0) = 0.01. The wave function, solution of the reduced 
NCWDW Eq. (EQ), is well-defined for the chosen values fl. 

Notice that the minimum of this potential is semiclassically stable. Following Ref. 
the probability of decay is given by c~ Se where Se is the euclidianized action for 
the bounce, i.e. for the motion of a particle travelling from xq = x(t ) (the position 
where the potential attains its local minimum) to x\ = x(t{) (the position where it attains 
V(xi) = after its local maximum) and arriving there with zero velocity. After lifting the 
potential so that V(xn ) — (which does not affect the probability of decay) one finds that 
Se = Ixd \j2V(x) dx 15]. Since during the relevant part of the motion the potential is 
essentially given by the second term in Eq. fl25l) . one can use values of Fig. [H which are 
fairly typical, to obtain Se — 9.5. 

In order to compute the partition function using the Feynman-Hibbs procedure one 
resorts to the eigenvalue problem of the harmonic oscillator. Thus we expand the exponen- 



ts 



tial term in the potential Eq. (1251) to second order in powers of the x — Xq variable. The 
minimum xq is obtained by solving the following equation 

dV, 



96^ exp (-2^0) - 277(7/20 - c) = . (26) 



dx 

The minimum is thus defined by the relation 

exp (-2V3x ) = (D- ^=x , (27) 

where £ = t//4a/3 and D = c/12. Moreover, one should impose that the second derivative of 
the potential function is positive in order to obtain the following condition to the minimum, 

6exp(-2v / 3xo)-C 2 >0^x <--^ln^^ . (28) 

We then get for the potential function V(x) to second order in x — x : 

V{x) = 48(6e- 2 ^ - ( 2 )(x - x ) 2 + A^e' 2 ^ - {rjx - c) 2 . (29) 

So, one can rewrite the ordinary differential equation resulting from the NCWDW equation 

as 

" l^S + 24 ( 6e ~ 2V ^° - C 2 )^ - x o) 2( P + [24e- 2 ^ - ^{r,x - c) 2 ]4> = . (30) 

Comparing Eq. (130]) with the Schrodinger equation of the harmonic oscillator we can 
clearly identify the potential and the energy of the system. Thus, the noncommutative 
potential required to compute the partition function is 

V NC (y) = 24(6e~ 2 ^ - CV , (31) 

where we have defined a new variable y = x — xq. The Feynman-Hibbs procedure allows for 
introducing quantum corrections to the partition function through the potential, i.e. one 
can admit a quantum correction for the potential given by 4J, 

i^BHrr/, 



24 



-V£ c (y) = 2f3 BH (Qe- 2 ^ - ( 2 ) , (32) 



where (3bh is the inverse of the BH temperature. The noncommutative potential with the 
corresponding quantum corrections then reads, 

U NC {y) = 24(6e- 2 ^ _ C 2) L> + Pbh\ (33) 



Finally, the noncommutative partition function is given by 

1 



'NC 



1 



cxp 



(34) 



\\ 48(6e- 2 ^° - C 2 ) Pbh 

We are now in conditions to compute the thermodynamic quantities of the Schwarzschild 
BH. Starting from the noncommutative internal energy of the BH, E NC = -^#-ln 

1 



d(3 E 



J NCi 



E 



NC 



P. 



BH 



+ 4(6e" 



-2-^/3x0 



BH 



(35) 



the equality E^c — M, allows for obtaining an expression for the BH temperature: 

3 — 

8(6e- 2 ^o _ (2) 

Inverting this quantity, one obtains the BH temperature. Dropping the term proportional 



l±(l-^(6e- 2 ^-C 2 ) Xl/2 



(36) 



to M 2 , as presumably M » 1, and considering the positive root, we obtain: 



T, 



BH 



c 2 



(37) 



Notice that this quantity is positive (cf. Eq. (|28|) ). To compare this result with the 
Hawking temperature Tbh = g^M' we should be cautious as the limit r\ — > is ill-defined. 
Indeed, our whole construction rests on the premise that we have a stable local minimum 
of the potential. However, this is only true provided rj > 0. We remark that our expression 
displays the same mass dependence as the Hawking temperature, and we write it as 

6(0 



T, 



BH 



M 



where h is a ( dependent constant: 

b(C) = 4(Qe- 2V * X0 - C 2 ) • 



(3? 



(39) 



Moreover, a simple numerical calculation reveals that we may recover the Hawking tem- 
perature even in the presence of the momentum noncommutativity for a specific value of 
rj. Indeed, equating Eq. (1371) with the Hawking temperature and using the stationarity 
condition (1270 . we obtain for c = 12D = 5.68 



x = 1.8478 r] = 0.025 . 



(40) 



And thus, since rj cannot be exactly equal to zero, we can regard 770 = 0.025 as a ref- 
erence value which yields the Hawking temperature, and as r\ increases we get a gradual 
noncommutative deformation of the Hawking temperature. 
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The entropy is calculated using the definition, S^c — ^Z NC + PbhEnc- Thus, the 
entropy for the phase-space noncommutative BH is given by: 



SBH = hi ^m + m)[ 1+ r- ^ 

M 2 




Neglecting as before terms proportional to M~ 2 , as M » 1, we finally obtain 

M 2 \/K0 

For the reference value r\ = r] , we have 6(Co) = We then recover the Hawking entropy, 
but also some "stringy" corrections 

Sbh = 4vrM 2 + In - ln(87rM) . (43) 

In summary, as in the case of the temperature, we obtain that the noncommutative entropy 
of the BH is the Hawking entropy plus additional contributions and some noncommutative 
corrections. 

As can be clearly seen, the relevant thermodynamic quantities are corrections of the 
Schwarzschild BH ones that depend explicitly on the momentum noncommutative param- 
eter 7] through 6(C)- If V = Vo, the Hawking temperature and entropy for the black hole 
are recovered. On the contrary, by inspection, if rj = 0, the potential Eq. (123]) reduces to 
a monotonous exponential term with no local minima. It is only if rj ^ that the poten- 
tial function acquires a quadratic term. This in turn ensures that the potential function 
has a local minimum and a local maximum. It is around the minimum of the potential 
that the calculations of quantum corrections can be carried out. Thus, the Feynman-Hibbs 
procedure is meaningful only in the presence of the momentum noncommutativity. 

IV. SINGULARITY 

In this section, we use Eq. (JTTJ) to extract some information about the essential r = t = 
singularity. It should be noticed that in order to employ the KS metric to describe the 
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interior of the Schwarzschild black hole the identification has to be imposed. Thus, we 
can see that for t = 0, Q — > +00 and (3 — ► +00. We are thus interested in studying the 
limit 

lim ^ c ,/? c ) , (44) 
\i c ,(j c — ++00 

where T/>(f) c , /3 C ) is a generic solution of Eq. (TP7|) . We seek for a representation of solutions 
of Eq. (117)) in terms of the eigenstates of A Eq. (fTB)) 

(45) 



V(a,&)= daC(a)M^cPc) 



where C(a) £ C and %l) a {Q. c , /3 C ) is of the form 



'a 



ip a (Q c ,l3 c ) = (p a I — I exp 



(46) 



and (j) a (z), z — ^c/a* satisfies Eq. ( I2ip . In the limit Q — ► 00, if we keep only the dominant 
terms, Eq. (j2~I)) reads: 

(47) 



4>'L( Z ) + fa* - a) 2 0a(^) = 
This equation can be rewritten for rj 7^ as 

r ^ 



77V o (z) = . 



(4* 



where the change of variables has been performed z = z — - and 4> a (x 
Eq. (T4H1) is analogous to the eigenvalue equation for an inverted harmonic oscillator. This 
Hamiltonian is self-adjoint in L 2 (IR 2 ), its spectrum is continuous and its zero eigenfunction 
(the solution of Eq. (|48)) ) displays the asymptotic form [if]] (for 77 7^ 0), 

1 



n/2 



exp 



±z^ 2 
2 



(49) 



and so, for all a, 
lim 



lim 



a, 

[z-- =0 

7/ 



lim ^ a (fi c ,/5 c ) = 

Oc,Pc— > + 0O 



(50) 



Therefore, it seems reasonable to expect that, for a suitable, although fairly general choice 
of the coefficients C(a) 

lim V(fio&) = 0, (51) 

S2 c ,Pc-*+oo 

which is a necessary condition to provide a quantum regularization of the classical singu- 
larity of the Schwarzschild BH. However, one should be cautious before running into the 
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conclusion that the probability of finding the BH at the singularity is zero. First of all 
the calculation of probabilities for general covariant systems is a subtle issue. In our case, 
given that the wave function is oscillatory in f3 C) it seems natural to fix a /3 c -hyper surf ace, 
corresponding to the introduction of the measure 8{(3 — /3 c )dj3dQ c in the probability distri- 
butions. The probability P(r = 0, t = 0) of finding the BH at the singularity would then 
be given by the expression: 

r+oo r+oo CV 

P(r = ,t = 0)= lim / |^r c ,/^|W ~ hm / \M-)\ 2 dn' c (52) 

n c ,/3 c ^+oo Jq c n c -^+ooJn c /i 

which, unfortunately, is divergent. This follows from the conclusion (which can be de- 
rived from the asymptotic expression) that the inverted harmonic oscillator displays non- 
normalizable eigenstates [if]]. Hence, the noncommutativity of the form ffl3l) cannot be 
regarded as the final answer for the singularity problem of the Schwarzschild BH. 

However, it is quite interesting to realize that an Hamiltonian H with the asymptotic 
form 

H = -— + V{x) , V(x) ~ -r] 2 x 2+2e (53) 
for some e > 0, displays zero energy eigenstates of the kind 17] 

V>(*) ~ ^7/2 ex P f^sk 5 **! , (54) 



which are normalizable. One may then speculate that a suitable perturbation of the non- 
commutative structure fflBl may lead to a NCWDW equation associated to a potential of 
the form V(x) ~ —rfx 2<yl+ ^ for some arbitrarily small e > 0. The solutions of this new 
NCWDW equation would then display zero probability (in the sense of Eq. (152]) ) at the 



singularity, thus solving this problem for the Schwarzschild BH 



1* 



V. CONCLUSIONS 



In this work, the NCWDW equation for a KS cosmological model was used to study 
the thermodynamic properties of the Schwarzschild BH and to obtain information about 
the structure of the singularity at t — r — 0. Working out the NCWDW equation we 
implemented the Feynman-Hibbs procedure identifying the quadratic part of the minisu- 
perspace potential. This allows for obtaining the partition function and hence the BH 
thermodynamic quantities. As we have shown, this is possible for a noncommutative BH 
only if noncommutativity in the momentum space is introduced. 

13 



The essential singularity r = t = is studied with the help of the NCWDW equation. 
In the limit of t = 0, or Q, (3 — > +oo, one arrives at the Schrodinger problem of the inverted 
harmonic oscillator, whose wave function vanishes at the limit Q,j3 — > +oo, but is not 
square integrable. Therefore, it is not possible conclude that the singularity problem is 
solved by the considered momentum noncommutativity. Nevertheless, we can foresee that 
a noncanonical form of phase-space noncommutativity may allow for a square integrable 
wave function. This possibility will be analysed elsewhere 18 1. 
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